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Abstract 

We consider a SDE with a smooth multiplicative non-degenerate noise and a pos- 
sibly unbounded Holder continuous drift term. We prove existence of a global flow of 
diffeomorphisms by means of a special transformation of the drift of Ito-Tanaka type. 
The proof requires non-standard elliptic estimates in Holder spaces. As an applica- 
tion of the stochastic flow, we obtain a Bismut-Elworthy-Li type formula for the first 
derivatives of the associated diffusion semigroup. 

1 Introduction 

In this paper we study the existence of a global stochastic flow of diffeomorphisms for the 
following stochastic differential equation in M. d 

k 

dX? = b(X?)dt + J2°i(X?)dWi, t>0, X$ = x, (1) 

i=l 

where Wt = (W - /, W t ) is a standard Brownian motion in M. k . We assume that the 
diffusion coefficients di : M rf — ► M d , i = 1,.. . ,k, are smooth and non-degenerate and we 
allow the drift term b : M. d — > M. d to be unbounded and Holder continuous. 

Following a common language, we say that equation (pQ) is weakly complete if there 
exists a unique global strong solution for every x E M. d , and that it is strongly complete if 
there exists a global stochastic flow of homeomorphisms. If the coefficients b and o~i are 
globally Lipschitz, then one has strong completeness (see |19] and |20j). 

Weak completeness is true under much weaker assumptions: for instance, when the 
coefficients b and o~i are locally Lipschitz continuous and have at most linear growth. In 
dimension one, these assumptions also imply strong completeness (see [19] and [20] ) but 
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in dimension larger than one there arc counterexamples, from |23j, even in the case of 
smooth bounded coefficients. These examples indicate that some form of global control at 
infinity on the increments of the coefficients is necessary. For (at least) locally Lipschitz 
coefficients, there are indeed positive results of strong completeness (see [TJ, [22], [21]). 

Strong completeness for non-locally Lipschitz coefficients can be established replacing 
the global Lipschitz condition on the coefficients with global log-Lipschitz type conditions 
(see [29], [33] . [9], [8]). Such log-Lipschitz conditions are stronger than the Holder conti- 
nuity. 

Many papers prove weak completeness for SDEs with non-locally Lipschitz continuous 
coefficients assuming a non-degenerate diffusion matrix a. First papers in this direction 
were [35] and [31] in which the method of the so called Zvonkin's transformation was 
introduced. More recent papers dealing with such approach are [13], [T7J, [32], [33] (see also 
the references therein). In the case of non-degenerate additive noise and time dependent 
drift b, the most advanced result (but see also the 1-dimensional results reported in [50] ) 
is [17] ; in such paper it is shown that it is sufficient to assume that b G L q (0,T; Lf oc (U. d ) ) 
with ! + |<l,p>2 and q > 2, plus a non-explosion condition, to get weak completeness. 
This result has been generalized in [32] to cover also the case in which a is variable, time- 
dependent and non-degenerate. We do not know about strong completeness under such 
weak assumptions. 

The contribution of the present paper is to prove strong completeness for SDEs with 
"locally uniformly ^-Holder continuous" drift b, for some 9 G (0, 1) (see ([3])), removing 
boundedness of b or additional regularity assumed in previous works. Also, we allow non- 
degenerate, bounded and C^(K. d , M d )-diffusion coefficients (<Ti)i = x,...,k- We point out that 
our result seems to be new even in the case of constant and non-degenerate (ci)i=i ... fe- 

In spite of the fact that b is not even differentiable, under the previous assumptions, we 
construct a stochastic flow of C 1 -diffeomorphisms (see Theorem [7J using the approach of 
[10] rather than the Zvonkin's transformation method used in the above mentioned works 
on strong completeness (we compare the two methods in Section 3). 

In [10] in order to study a linear stochastic transport equation with a bounded vector 
field b(t, x) which is Holder continuous in x, uniformly in time, we have showed that if in 
dH) a = (<Tj) is constant and non-degenerate and b = b, then there exists a stochastic flow of 
C 1 -diffeomorphisms. This result can be extended without difficulties to the case in which 
a is not constant, bounded, non-degenerate, and time-dependent (see |34j where this case 
is investigated by the Zvonkin's transformation or Remark [9] where we show such result 
following the approach of |10j). 

In the present situation, since our b is unbounded, we need new global regularity results 
in Holder spaces for the solution u of the elliptic equation 

\u(x) - ^Ti{a(x)D 2 u(x)) - b(x) ■ Du{x) = b{x), x G R d , (2) 

to be interpreted componentwise, where A > is large enough, a(x) = a(x)a*(x) (cr*(x) 
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denotes the adjoint matrix of cr(x)). The study of this equation will be the subject of 
Section 2 of the present paper. The required estimates are not covered by recent papers 
dealing with elliptic and parabolic equations with unbounded coefficients (compare with 
@]) [I]) PS] and the references therein). To obtain such result we prove a crucial Lemma 
[3] concerning estimates on the derivatives of the associated diffusion semigroup when it is 
applied to unbounded functions /; in its proof we also use an argument from the proof of 
\27\ Theorem 3.3]. In Remark [10] we show a possible extension of our Theorem [7] to the 
case in which b and a are time-dependent. 

We finish the paper by showing that a Bismut-Elworthy-Li formula holds for the diffu- 
sion semigroup associated to ([I]) (see Theorem II ip . Under the poor regularity of b assumed 
here, this result is new. Bismut-Elworthy-Li formula requires a suitable form of differen- 
tiability of the solution of ([1]) with respect to the initial condition x; we have this result as 
a byproduct of our Theorem [7] on existence of a differentiable stochastic flow. 

Notations and assumption The euclidean norm in any k > 1, will be denoted by 
| • | and its inner product by • or (•, •). For 9 £ (0, 1), we define the set C (R ; k, d > 1, 
as set of all vector-fields / : M. d — > R fc for which 

[f\e := sup — j - 9 — < oo. (3) 

,\x-y\<l \ x V\ 



These are the "locally uniformly 9-Hdlder continuous" vector fields mentioned in the intro- 
duction. The function / (x) = \x\ e is a classical example. We let 

r/i \f( x )-f(y)\ / (A \ 

[f]0 ' i:= SUP , flr ^Vlr ? ;h <0 °' (4) 

x^ y eR d \\ x y\ v \ x y\) 

where a V b = max(a, b), for a, b £ M. By a simple argument we have [f]g < [f]g i < 2[/]g, 
so in particular functions in C (R d ;R ) have at most linear growth. The set C e (R rf ;R fc ) 
becomes a Banach space with respect to the norm 

11/11* = || (1 + 1 ■ I) -1 /Olio +[/]«• 

where || • ||o denotes the supremum norm over R rf . We say that / G C n+e (]R rf ; R fc ), n > 1, if 
/ G C e (M. d ; R fc ) and moreover, for alH = 1, .. . ,n, the Frechet derivatives D i f are bounded 
and ^-Holder continuous. Define the corresponding norm as 

n 

+ Y f \\B l f\\ o + [D»f] . (5) 

i=l 



If R k = R, we simply write C n+e (R d ) instead of C n+9 (R d ;R), n > 0. C£ +e (R d ; R fc ) is 
the subspace of C n+e (R d ; R fc ), consisting of all bounded functions of C n+e (R d ;R k ). In 
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particular, C®(M. d ) is the usual Banach space of all real bounded and ^-Holder continuous 
functions on R d (cf. [Ml). C^(R ;R ) is the space of all bounded functions from R d into 
R fc having also bounded derivatives up to the order n > 1 and we set C£(R d ; R) = C£(M. d ). 
Finally, we say that / : R rf — > R d is of class C n,a , n > 1, a G (0, 1), if / is continuous on 
R rf , n-times differentiable and the derivatives up to the order n are a-H61der continuous 
on each compact set of M. d . 

Throughout the paper we will assume a fixed stochastic basis with a d-dimensional 
Brownian motion (SI, (Ft) , J 7 , P, (Wt)) to be given. Denote by T s ,t the completed <r-algebra 
generated by W u — W r , s < r < u < t, for each < s < t. 

On equation (TTJ), we will consider the following assumptions. 

Hypothesis 1 There exists 9 G (0, 1) such that b G C 9 (R d ;]R d ). 

Hypothesis 2 The diffusions coefficients o~i : M. d — > W d , i = 1, . . . , k, are bounded func- 
tions of class C^{R d ,R d ). 

Hypothesis 3 Consider the d x k matrix o~(x) = (o~i(x)), and its adjoint matrix cr*(x), 
x G M. d ; we assume that, for any x G M. d , there exists the inverse of a(x) = o~(x)a* (x) and 

1 1 a- 1 1 1 o = sup ||a _1 (x)|| < oo (6) 

x£R d 

(||a _1 (x)|| denotes the Hilbert-Schmidt norm of the d x d symmetric matrix a~ 1 (x)). 

2 Regularity results for the associated elliptic problem 

2.1 Estimates on the derivatives of the diffusion semigroup 

Here, we consider the SDE ([1]), assuming that a satisfies Hypotheses [2] and [3] and imposing 
in addition that 

b G C 3 (R d ;M d ) with all bounded derivatives up to the third order. (7) 

Clearly this is stronger than Hypothesis [T] but b is not assumed to be bounded. 

Let (Pt) be the corresponding diffusion semigroup, i.e., for any g : M. d — > R Borel and 
bounded, 

P t g(x)=E[g(Xf)], x G R d , t > 0, 

where (Xf) is the unique strong solution to ([I]) under (J7]). 

In our next result, we will prove estimates on the spatial derivatives of Ptf, t > 0, 
assuming that / G C e (R d ). To this purpose, we will use the so-called Bismut-Elworthy-Li 
formula (see (fT2|) ) for the spatial derivatives of Ptf (cf. [6]). 
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Let us comment on such formula. Probabilistic formulae for the spatial derivatives of 
Markov semigroups have been much studied for different classes of degenerate and non- 
degenerate diffusion processes even with jumps (see [3], [21], [6], [5], [3], [12] [26], [27] . 
|34j and the references therein). The martingale approach of [6] mainly works for non- 
degenerate semigroups (but see also [H Chapter 3] and [31]); it has been also used for 
some infinite dimensional diffusion processes (see [5] and [4]). On the other hand, in 
case of degenerate diffusion semigroups, more complicate formulae for the derivatives can 
be established by Malliavin Calculus (see [3], [21], [12] and [26]). Some applications to 
Mathematical Finance are given in |14j . 

The next lemma is of independent interest since the function / in (|8|) is not assumed 
to be bounded (compare with [U Chapter 1] and [I] Chapter 6]). 

Lemma 4 Assume Hypotheses^ and\^ and condition (|7|). There exist constants Cj > 0, 
Mj > 0, j = 1,2,3 (cj and Mj depends on 6, \\a~ l \\o, d, ||a||o and on the supremum norms 
of derivatives of a and b up to the order j), such that, for any f G C e (R d ), t > 0, it holds 

c ■ t 

llDtPtfh^MjWe^p, t>0, for j = 1,2,3. (8) 

Proof. / Step. First note that E[sup te[0]T ] \Xf\i] < C T (1 + \x\i), for any T > 0, x G R d , 
q > 1 (see, for instance, [191 Chapter II]). 

It is also known that, for any t > 0, the mapping: 

x i — y is three times Frechet differentiable from M d into L 2 (Q) (9) 

(see [4] Section 1.3] which contains a more general result). Let us write the Frechet deriva- 
tives: 

th(x, h) = D x (Xf)[h], e t (x, h, k) = D 2 x (Xf)[h, k], -fo{x, h, k, I) = D 3 x (Xf)[h, k, l], 

for any x, h,k,l G M. d . These derivatives satisfy suitable stochastic variation equations (see 
[191 Chapter II]). We only write down the variation equation for ijt = ijt(x, h): 

d Vt = Db{Xf) m + Da{Xf) m dWu Vo = h. 

Using standard estimates, based on the Burkholder inequality, we get that, for any p > 1, 
that there exist positive constants C and c (depending on p, ||-Dfr||o an d ||-Dc||o) such that, 
for any x G R d , h G R d , 

E\m(x,h)\ p < C\h\ p e ct , t > 0. (10) 
In a similar way, using the second and third variation equations, we obtain the estimates: 

m(x,h,k)\? <c 2 \hnk\^\ (ii) 
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E|V>tO, h, k, l)\ p < C 3 \h\ p \k\ p \l\ p e fi3 ', t > 0, 

for any x, h,k,l E (with positive constants C\ and c\ which depend on p and on the 
supremum norms of the derivatives of b and a up to the order i, i = 2, 3). 

// S^ep. Arguing similarly to [4j Section 1.5] one can prove that, for any / E C (R d ), 
t > 0, the map: x h- > P t f(x) is differentiable on ]R rf and, moreover, we have the following 
Bismut-Elworthy-Li formula: 

(DP t f(x),h) =K\f(Xf)J 1 (t,x,h)], x, h E R d , t > 0, where (12) 

J\t,x,h) = \ [\a*(X*)a- 1 (X*) Vs (x,h),dW s ). 
z Jo 

Note that formula (|12|) is first proved for bounded / E C|(M rf ). Then a straightforward ap- 
proximation argument shows that (|12j) holds even for (a possibly unbounded) / E C (R d ). 
However, to be precise, in [3], it is assumed that <r(x) is an invertible d x d matrix and so 
the expression of J 1 in [4j Section 1.5] contains a (X®) instead of our a*(Xf) a^ 1 (X^). 
We briefly explain why (|12p holds following the proof of |25|, Theorem 5.1]. We only discuss 
the crucial point of the argument which is needed to get (|12p when / E C|(]R rf ). One has 
by the Ito formula 

f(Xf) = P t f(x) + [ (DP t _ s f{X x s ),a{X x s )dW s ). 
Jo 

Multiplying both terms of the identity by the martingale 

K t = f\a*(X*)a- 1 (X*)r ]s (x,h),dW s ), 
Jo 

and taking the expectation, one arrives at 

E[f(Xf)K t ] = (\[{DP t ^ s f{X x s ), Vs {x : h))\ds = t(DP t f(x),h). 
Jo 

Thus (|12|) is proved. 

Now the problem is to show that, for f E C d (M. d ), t > 0, the map: x i— > (DPtf(x), h) is 
a bounded function (we cannot use as in the boundedness of f). 

By using (fT0|) . we get easily that there exist C\ > depending on ||a _1 ||o, ||-D&||o an d 
||Dcr||o such that 

p Cit 

E\J l {t,x,h)\ 2 < -i^|/i| 2 , t>0. (13) 

Now we prove the crucial estimate of the first derivative in ([8]) . We use an argument from 
the proof of \27\ Theorem 3.3]. Introduce the deterministic process 

't 



Y t x = x + I b{Y s x )ds, t > 0, x E R d , 
Jo 
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which solves Y t x = b(Yf), Y X = x. Using that a is bounded and applying the Gronwall 
lemma, we find, for any q > 1, 



E\Xf - Y t x \ q < Mt q/2 e Clt , t > 0, x G R d , 
where M depends on \\<j\\o and q and c\ on ||-D6||o and q. Since 

E[f(Yf)J l (t,x,h)] = f(Y?)(D(P t l)(x),h) = 0, t >0, /i G M d , x G 
we have (see also fl3J)) 

|<£>P t /(s),Zi>| = |E[(/(Xf)-/(F^))J 1 (t,x,/ i )] 

<2\f] e E[(\Xf-Y t *\°V\X?-Y t *\) \JHt,x,h)\] 
< 2[f] e (E[\Xf - Y*\ 20 V \Xf - Yf\ 2 }) 1/2 (nj\t,x,h)\ 2 ) 112 , 



(14) 



(15) 



t > 0. Using that a V 6 < a + b, a, b > 0, and the previous estimates (|1U|) and (|14|) . we find 

J't C , e c't 



\h\, t > o, xer, (16) 



|(DP«/(x),fc>| < C"[/] e (^ 2 + t 1 /2 ) _| /l | < [f]e ^ 

where C and c depend on ||cr||o, 1 1 a. 1 1 1 o 5 ||-^ cr ||oj H^llo an d 0. 

Let us consider the remaining estimates in ([8]). We have, using the semigroup law, 
P t f = P t/2 {P t/2 f) and so (cf. U formula (1.5.2)]), for any x, h, k G R d , t > 0, 

(D 2 (P t f)(x)k, h) = DM{P t/2 f){X { t ] 2 ) j\t/2, (•), h)] ) (s) 



E 



(DP t/2 f(Xf /2 ), rh /2 (x, k)) J\t/2, x, hj\ + E[p t/2 /(Xf /2 ) D k J x (tj2, x, h) 

= r 1 (t,x) + r 2 (t,x), 



where denotes the directional derivative along the vector k (indeed, for any fixed t > 
and h G R d , the mapping: x \— ► J l (t/2,x,h) is Prechet differentiable from K rf into L 2 (J7); 
this follows easily, using ©, CHD, (JTOJ) and ([IT])). We have 

D k J 1 (t/2,x,h) = - / (^(^^[^(x^^a- 1 ^)^^,^,^) 
t ./n 



2 f t/2 

7 



<<t*(X*) a" 1 ^) DflK)^^ k)} a~\X x s ) Vs (x, h),dW s ) 
c Jo 
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Using the Schwarz inequality, (fT3|) and 

/-ill „c'Ot 

s ^{n{DP t/2 f{x't /2 )^ t/2 {x,k))\ 2 ) l/2 < ifhjrj^lkl 



we get immediately |Ti (t, x) | < M[f]g t (2-e)/2 \ h\\k\, t > 0, x £ M d . To estimate T 2 , first 
note that, by taking / = 1, 

= {D\P t l)(x)k,h) = + E[D k J 1 (t/2,x,h)], 

for any x, h,ke R d . We find (ar guing similarly to (I15p ) 

T 2 (t,x) =E[(P t/2 f(X? /2 ) - P t/2 f(Yf /2 )) D k j\t/2,x,h) . 

Since 

\P s f(x) - P s f(y)\ < E\f(Xf) - f(Xy)\ < 2\f] e E[(\Xf - Xjf + \X* a - X^) 

< 2[f] e M(\x-y\ 9 + \x-y\)e c ' s , s>0, x, y G R d , 
we find, for any x G R rf , t > 0, 

|r 2 (t,x)| <2Me ci */ 2 [/],E[(|xf /2 -y^ 2 | e +|^ /2 -y^ 2 |) ia^A^)!]. 

< 2Me^ [/]e ( E [|Xf /2 - Y^ 2 | 2e + |Xf /2 - ^ 2 | 2 ]) 1/2 (E|X) fc J 1 ^, s, Ml 2 ) 1 / 2 

where Ci and c\ depends on ||<r||o, ||a _1 ||o, ||-Dcx||o, ||-D 2 cr||o ||D6||o, ||D 2 6||o and 6. We have 
so obtained estimate in (J8|) corresponding to j = 2. 
The estimate for j = 3 follows in a similar way. ■ 



2.2 The main regularity result 

With respect to the previous section, here we consider the elliptic operator 
Lu(x) = ^Tr(a(x)D 2 u(x)) + b(x) ■ Du(x), x G R d , 

with a{x) = a(x)a*(x), assuming Hypotheses [Q [2] and O 

The next result provides new estimates for L in Holder spaces. These estimates are 
not covered by recent papers dealing with elliptic and parabolic equations with unbounded 
coefficients, due to the fact that in our case also / can be unbounded (compare with [1], [I], 
|18j and the references therein). 
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Theorem 5 Let E (0,1). For any 0' E (0,9), there exists Ao > (depending on 0,0', d, 
[b] g , \\o-\\o, Hct" 1 ||o, ||- Dfe cr||o, k = 1,2,3) such that, for A > A , for any f E C e (R d ), the 
equation 

Xu-Lu = f (17) 
admits a unique classical solution u = u\ E C 2+e ' (R d ) for which 

\\u\\ 2 +o> = IK) (1 + I • I) -1 No + 11-DttHo + \\D 2 u\\ + [D 2 u] e i < C(X)\\f\\ e (18) 
with C(A) (independent on u and f) such that C(A) — > as A — > +oo. 

Proof. Uniqueness can be proved by the following argument (cf. [16| page 606]). Consider 
7](x) = a/1 + \x\ 2 , x £ R d . 

Defining u = vrj, we obtain an elliptic equation for the bounded function v, i.e., 

Xv(x) - -Tr(a(x)D 2 v(x)) - (b(x) + a ^ D ^ ) . Dv ( x ) 
2 r]{x) 

17V(.(,)i>V«)) + D^Wx /w_ 



r/(x) ?y(x) 7 7?(rr) ' 

x S Note that v has first and second bounded derivatives. For A large enough (de- 
pending on ||o"||o and ||^||o), uniqueness of v follows by the classical maximum principle. 

Now we divide the rest of the proof in some steps. 

Step I. We assume in addition that b E C 3 (M rf ,]R ci ) and has all bounded derivatives up 
to the third order (but it is not necessarily bounded). We prove that, for sufficiently large 
A > 0, there exists a unique solution u = u\ E C 2+ (R d ) to the equation 

Xu-Lu = f E C e (R d ). 

Moreover there exists C (independent on u and /) such that 

hh+e < C\\f\\ e . (20) 

Estimates (|20p are new Schauder estimates since / is not assumed to be bounded (compare 
with @] and [I]) 

We consider the function 

u(x) = / e' xt E[f(Xf )}dt = / e~ xt P t f(x)dt, x E R d , (21) 
Jo Jo 

where (Xf ) is the solution of (pQ) and show that, for A large enough, u is a C 2+6, (M d )-solution 
to our PDE. 
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Using that E\Xf - Xf\ < Ce ct \x - y\, t > 0, x, y E R d , we find 

\u(x) - u{y)\ < c[f] e ,i (\x - yf V \x - y\), x, y G R d , 

and also ||u(-) (1 + | • |) _1 ||o < C||/(") (1 + I " l) _1 ||cb f° r A large enough. 
By Lemma S] we get, for A large enough, 

\\Du\\ +\\D 2 u\\ < C[f]g. 

To estimate the second derivatives of u, we proceed as in \27\ Theorem 4.2]. We have, for 
any x, y £ M. d with \x — y\ < 1, 

\D 2 u{x) - D 2 u{y)\ = Vl e- xt \D 2 P t f(x)-D 2 P t f(y)\dt 
Jo 

POO 

+ e~ xt \D 2 P t f(x)-D 2 P t f(y)\dt 

J\x-y\ 2 

< c"\x - y\ e [f} e + C\x - y\ [f] [°° e~ xt dt < c'[f] e \x - yf. 

J\x—y\ 2 " 

It remains to check that u is a solution. This is not difficult thanks to Lemma [4] (see, for 
instance, [H Chapter 1] or argue as in [27\ Theorem 4.1]). 

Step II. Under the assumptions of Step I, for any a £ (0, 9), we have 

H 2+ a<C(A)||/|k (22) 

with C(A) — > 0, as A — ► +oo. This is clear if we replace ||u||2+ a with \\u(-) (1 + | • |) _1 ||o 
+||-Dit||o + ||-D 2- u||o- Therefore, we only consider [D 2 u] a . 

Combining the interpolatory estimate: [v] a < C||f Ho"" ll-^-Mloi v ^ (where 
C = C(d), see [151 Section 3.2]) with estimates of Lemma H] corresponding to j = 2,3, we 
find, for any t > 0, 

„C4* 

[D 2 P t f] a < C\\D 2 P t f\\l~ a \\D 3 P t f\\% < C A [f] e — , 
with 7 = < i ( s i nce a <9). It follows 

r+oo (c 4 -A)t 

[D 2 u] a < C A [f] e / —p=—dt < C 5 [f] e (X - c^' 1 . 
Jo t ' 

The assertion is proved. 
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Step III. We require that b G C 9 (R d ,IR rf ) as in Hypothesis [T] and prove the following 
a-priori estimates: if A is large enough and u G C 2+e (K d ), < 9' < 9, is a solution to 
Xu-Lu = f e C (M. d ), then 

IK) (1 + I • ir'llo + \\Du\\ + \\D 2 u\\ e > < K(X)\\f\\e, (23) 

with K(X) — ► 0, as A — > +oo. 

To prove the estimate we introduce p G Co°(W i ), < p < 1, = p(—x), for any 
£ G J" (ix = 1. Moreover, b * p indicates b convoluted with p. 

Write Xu(x) - \Tr (a(x)D 2 u(x)) — (b * p){x) ■ Du(x) = f(x) + ({b — (b * p)) (x), Du{x)). 
It is easy to see that b* p (even if it can be unbounded) is a C°°— function with all bounded 
derivatives. Moreover, there exists C = C(9, Dp, D 2 p, D s p) > such that 

\\D k (b* P )\\ <C[b] e , k = 1,2,3. (24) 

The function b — (b* p) is bounded and we have 

\\b-(b*p)\\ <C[b} e . 

It follows that b - (b * p) G C%(R d , R d ). Applying Step II, we find that 

IH| 2+ e'<C(A)||/|| fl + C(A)||<6-(6*p), J D«)|| fl (25) 

with C(A) 0. Using that 

||(6- (b* p),Du)\\ e < c[b]e\\Du\\ + c[b] g \\Du\\ g < c[b} e \\u\\ 2+ e > , 

for some constant c depending on 9, we rewrite ([25D : 

\\u\\ 2+e > < C(X)\\f\\ e + C(X)clb} e \\u\\ 2+ e'. 
Choosing Ao > such that C(A) < ^r^, for A > Ao, we find, with u = u\ 

(1 - C(X)c[b} e )\\u\\ 2+e , < C(X)\\f\\ e . (26) 

Defining K(X) = i_c(\)c{b] g ' we the assertion. 

Step IV. We show that for A > Ao (see Step III) there exists a classical solution u = u\ G 
C 2+e (M. d ) to (|17p . This assertion will conclude the proof. 

We fix A > Ao- To prove the result, we will use the continuity method. To this purpose, 
using the test function p of Step III, we consider: 

Xu(x) - ^Tr(a(x)D 2 u{x)) -(!-<*)(&* p){x) ■ Du{x) - 5b{x) ■ Du{x) = f{x), (27) 
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x G R , where 5 G [0, 1] is a parameter. Let us define 

r = {5 G [0, 1] : eq. (P27]) has a unique solution u = u 5 G C 2+9 '(R d ), for any / G C e (R d )}. 

T is not empty since G T by Step I. Let us fix Jo G T and rewrite equation ([27]) corre- 
sponding to an arbitrary 5 G [0, 1] as 

Xu(x) - ^Tr(a{x)D 2 u(x)) - (1 - 5 )(b * p)(x) ■ Du(x) - 5 b{x) ■ Du(x) 

= f(x) + [S-S ] (b- b*p)(x)-Du(x). 

Introduce the operator T : C 2+e '(R d ) C 2+e '(R d ). For any v G C 2+e '(R d ), Tv = u is the 
(unique) C 2+e (R d )-function which solves 

\u{x) - ^Tr(a{x)D 2 u(x)) - (1 - 5 )(b * p)(x) ■ Du(x) - 5 b(x) ■ Du(x) 

= f{x) + [6 -So] (b- b*p)(x)-Dv(x). 
Using the a-priori estimates ()26|) . we get that 

\\Tv - Tw\\ 2+e > < 2K(\)\S-S \ [b] e \\v-w\\ 2+e >, v,w G C 2+e '(R d ). 

Choosing \5 — 5q\ small enough, the operator T becomes a contraction on C 2+e ' (R d ) and it 
has a unique fixed point which is the solution to (|27|) . Therefore for \5 — 5q\ small enough, 
we have that 5 G T. A compacteness argument shows that T = [0,1]. The assertion is 
proved. ■ 

3 Different iable stochastic flow 

Given x G R d , consider the stochastic differential equation in R d : 

dX t = b (X t ) dt + a(X t )dW t , X s = x, t>s>0. (28) 

As already mentioned our key result is the existence of a differentiable stochastic flow 
(x,s,t) i — > ip S) t(x) for equation (|28[) . Recall the relevant definition from H. Kunita [19j : 

Definition 6 A stochastic flow of diffeomorphisms (resp. of class C 1,a ) on the stochastic 
basis (fi, (.Ft) P, (W*)) associated to equation (28\) is a map (s,t,x,uj) t— > <f> s ,t(x) (u), 
defined for < s < t, x G R d , lo G Q with values in R d , such that 

(a) given any s > 0, x G R d , the process X s,x = {X^ x (u>) , t > s,u> G O) defined as 
X^' x = (j) St t(x) is a continuous F s ^-measurable solution of equation I128\) ; 
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(b) P-a.s., for all < s < t, (j) s t is a diffeomorphism and the functions (j) s t(x), 4> st (x), 
D(j) s j(x), D(f)~l(x) are continuous in (s,t,x) (resp. of class C l,a in x uniformly in 
(s, t), /or <s <t <T, with T > 0); 

(c) P-a.s., 4> s ,t(x) = (/) u ,t(<l)s,u{x)), for all < s < u < t, x G M. d , and 4> SjS (x) = x. 

Starting from the work of Zvonkin, an important approach to the analysis of SDEs 
with non-regular drift is based on the transformation ^ : M. d — > M. d , solution of the vector 
valued equation 

-^ + »( = 0on [0, T] , ^ T (x) = x 

where ^t(x) = ^f(t,x) and [0,T] is a time interval where the SDE is considered. At time 
T, the solution is an isomorphism by definition; one has to prove suitable regularity and 
invertibility of ^ t for t G [0, T}. Then Y t := $> t (X t ) satisfies 

dY t = DV t (y^ 1 (Y t )) a i^ 1 (Y t )) dW t . 

The irregular drift has been removed. This approach, although successful (see [2], [13], 
[17] . [32j . [34]), raises two delicate questions: i) one has to deal with unbounded initial 
conditions; ii) one has to prove some form of invertibility. 

We propose a variant, based on the same operator L but on the vector valued equation 

\ip - Lip = b 

(under other assumptions one can treat also the time-dependent case through the parabolic 
equation \ip t — — Lip t = b, see [HI]). We find it more tractable than the case of unbounded 
initial condition; and we translate the difficult invertibility issue in the smallness of the 
gradient of the solution, obtained by means of a large A. When the gradient of ip is less 
than one, the function ^ (x) = x + tp (x) is invertible and the process Yt := ^ (Xt) satisfies 

dY t = (m- 1 (Y t )) a (V- 1 (Y t )) dW t + Xip (^ (Yt)) dt. 

So, at the end, the transformed equation has the same degree of difficulty as in the case of 
the Zvonkin's transformation. 

Theorem 7 Assume Hypotheses\^\^\^and fix any 0" G (0,6). Then we have the following 
facts: 

(i) (pathwise uniqueness) For every s > 0, x G M. d , the stochastic equation A28\) has a 
unique continuous adapted solution X s,x = (Xp x (uj) ,t > s, to G $7). 

(ii) (differentiate flow) There exists a stochastic flow <j) = (<fi s ,t) °f diffeomorphisms for 
equation [28\). The flow is also of class C ' . 
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(Hi) (stability) Let (b n ) C C 9 (M rf ,]R d ) and let (4> n ) be the corresponding stochastic flows. 
Assume that there exists b G C e (R d ,R d ) such that b n - b £ C%(R d , R d ), n > 1, and 
\\b — b n \\ c e — > as n — > oo. If (j) is the flow associated to b, then, for any p > 1, 



T > 0, 



\4>^ u (x)-(t> su (x)\ p 
km sup sup E[ sup — ■ j- ' i rhz ] = 0. (29) 

io<s<t ue[s,T] U + Flr 



sup sup sup E[ sup \\D^ u (x)\\ p ] < oo, (30) 

ngN x£M. d 0<s<T «e [s,T] 

Jim sup sup E[ sup ||D<# u (x) - D^ jU (x)|| p ] = 0. (31) 

iQ<s<T ue[s,T] 



(\\ ■ || denotes the Hilbert- Schmidt norm). 

Proof. Step 1 (auxiliary elliptic systems). Let us choose 9' such that < 6" < 9' < 9. 
For a fixed A > Ao > (see Theorem [5]) we consider the unique classical solution 
ip = ip\ G C 2+e '(R d ,R d ) to the elliptic system 

XiPx - L^a = b, (32) 

where 

Lu(x) = -Tr{a{x)a* (x)D 2 u(x)) + b(x) ■ Du(x), 

for any smooth function u : R d — ► R d (clearly (|32|) has to be interpreted componentwise). 
Define 

$x(x) = x + Va(z)- 
Similarly to [10l Lemma 8] we have 

Lemma 8 For A large enough, such that \\D7pxW0 < 1 (see Theorem^), the following 
statements hold: 

(i) ^ \ has bounded first and spatial derivatives and moreover the second (Frechet) deriva- 
tive D 2 ^>x is globally 9' -Holder continuous. 

(ii) ^x is a C 2 -diffeomorphism ofR d . 

(Hi) 1 f A ~ 1 has bounded first and second derivatives and moreover 

D*- x \y) = Y,(- D ^x\y)))\ y^ d - (33) 

fc>0 

In the sequel we will use a yaiue of A for which Lemma [3 holds and simply write ip and 
VP for ipx and ^x- 
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Step 2 (conjugated SDE). Define 



b(y) = X^-^y)), a{y) = D*{*-\y))a{y-\v)) 
and consider, for every s > and y G R rf , the SDE 

y f = y + T a(y u )dW„ + [ b(Y u )du, t > s. (34) 

J s J s 

This equation is equivalent to equation (|28p . in the following sense. If Xt is a solution to 
(j28[) . then Yj = ^(X^) verifies equation (|34p with y = ^>{x): it is sufficient to apply ltd 
formula to *S?(Xt) and use equation (f32T) . 

Viceversa, given a solution Y t of equation (f34")) . let X t = \I r ~ 1 (Yt), then it is possible to 
prove by direct application of Ito formula that Xt is a solution of (|28|) with x = 
This is not very important since below we will obtain this fact indirectly. 

Step 3 (proof of (i) and (ii)). We have clearly b and a G C 1+e ' (with first order 
derivatives bounded and in C® ) so that, in particular, they are Lipschitz continuous. 
By classical results (see \19\ Chapter 2]) this implies existence and uniqueness of a 

strong solution Y of equation (|34p and even the existence of a C 1 ' stochastic flow of 
diffeomorphisms ip S) t associated to equation (p4l) . 

The uniqueness of Y implies the pathwise uniqueness of solutions of the original SDE (HD 
since two solutions X, X give rise to two processes Yt = ^(X t ) and Y t = ^(X t ) solving (f3"4"|h 
then Y = Y and then necessarily X = X. By the Yamada-Watanabe theorem pathwise 
uniqueness together with weak existence (which is a direct consequence of the Girsanov 
formula) gives the existence of the (unique) solution {Xf)t> s of eq. (pQ) starting from x at 
time s. Moreover setting 

4>s,t = I'" 1 ° <Ps,t ° * 
we realize that <p St t is the flow of (pQ) (in the sense that Xf = <fr s t(x), P-a.s.). 

Step 4. (proof of (iii)). Let ip n and ip be the solutions in C 2+e ' (R d ;U d ) respectively 
of the elliptic problem associated to b n and to b G C e (R. d ;W i ). Notice that we can make a 
choice of A independent of n. We write 

A (-0™ —ip) — L (ip n - ^) = (b n -b) + (b n - b) ■ Dip n , n > 1. 

By Theorem[5]we have sup n>1 H^nll^+e' < C < oo. Since b — b n is a bounded function, by 
the classical maximum principle (see [15]) we infer also that tp — ip n is a bounded function 
on R d and 

C + 1 

||V'-'0n||o< — — ||6-fen||o) n>l. (35) 
It follows that ^ - 4) n G C 6 2+0 '(R d ;R d ) and ||^ - ip n \\ C 2+e> -> as n -» oo. 
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Fix p > 1 and consider the flows 99™ t = ^ n o (j>™ t o (^ n ) 1 which satisfy 

<PlM = V+ f t b n o^(y)du+ fa n o^(y).dW U) (36) 



We have a n — > a and b n — ► b, as n — > 00, in C 1+e ' (R d ;R dxk ) and C 1+d '(R d ; R d ), respec- 
tively. By standard arguments, using the Gronwall lemma, the Doob inequality and the 
Burkholder inequality (compare, for instance, with the proof of |191 Theorem II. 2.1]) we 
obtain the analog of (f29j) for the auxiliary flows ip™ t and ip S) t\ 

\V n su {x) - ip s ,u{x)\ P 

lim sup sup E\ sup — '■ — : — ^ ] = 0. (37) 

n - >oo xeR«'0< a <T ue[s,T] { 1 + \x\) p 

We can also prove the inequality 

sup sup sup E[ sup \\D<p™ u (x)\\ p ] < 00, (38) 

nGN x eK d 0<s<T ue[s,T] 

for D(p™ t (y), using the fact that the stochastic equation for Dip™ t (y) has the identity as 
initial condition and random coefficients Db n (0™ u ) and Da n which are uniformly 

bounded functions (since ||-D6 n ||o +||-D<j n ||o < C, uniformly in n). 

To prove (|3U|) is then enough to estimate D4> n su using (f55j) . the uniform boundedness of 
the derivatives of ^ n and its inverse (note that the uniform boundedness of the D(^ n )~ 1 
can be proved by (f33j) ). 

To prove (f29j) we remark that to estimate the difference ip™ t (ty n (x)) — (p St t(^(x)) we 
can split it as (^™ t (^ n (x)) — (p Sjt (^ n (x)) + ip S)t (^ n (x)) — (p Stt (^(x)). The two differences can 
then be controlled by 

E[ sup |<tf u (tt"(x)) - ip s , u (V n (x))\ p ] < «n (1 + |* n (x)|) p < a n (1 + \x\f, 

s<u<T 

(where a n = sup xeMd sup < s < T E[sup ue[sT] ^"'^^y^ ] and lim^oo a n = 0) and by 
E[ sup \v s>u (* n (x)) - <p s , u (*(x))\ p ] < su P E[ sup \\D^ u (z)\\ p ]\* n (x)-*(x)\ p 

s<u<T z t=Rd s<u<T 

< c\\y n -^\\ p , 

with lim^oo \\V n - ||o = lim^oo ||^ n - Vllo = (see (05}). 

Finally, one has to check that (*I' n )~ 1 converges to \E'~ 1 in the supremum norm. This 
follows from the inequality 

sup |(¥ B ) -1 (i/) - ^-^y)! < sup |(* n )- 1 (* n (x)) - ^-^^(x))! 
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< sup - < ||£># _1 ||o \\^ - * n || , 



which tends to 0, as n — > oo. 

Arguing as in the proof of |19^ Theorem II.3.1], we get the following linear equation for 
the derivative Dcj) s j(x) 

[D^{^ t {x))]D^ t (x) = D^(x) + f\DH{^ u {x))]D^ u (x)a( ( j )S)U (x))dW u 

t t ( 39 ) 

+ / D*(^ S)U (x))[D ( 7(0 s , u (x))] J D0 S)U (x)dTy u -A / [D$(4> at *(x))]D<t> a>u (x)du, 



< s < t < T, x G R d . From the fact that lim^oo \\tp n - ^|L 2+ s' = together with (130]) 
and (|39|) . we finally obtain 

lim sup sup E[ sup \\D^ u (x)-D<p S)U (x)\\ p ] = 0, (40) 

n ^°°xGlR d 0<s<T ue[s,T] 

which concludes the proof. ■ 

We consider now two possible extensions of Theorem [7J to the case when coefficients b 
and <7j are time-dependent continuous functions defined on [0,T] x M. d , i.e., we are dealing 
with 

k 

dXf = b(t,X?)dt + J2 (7 i( t i X t) d Wi, te[0,T], X = x. (41) 

i=l 

Remark 9 Let us treat the case in which also b is bounded. Following [10] . an analogous 
of our Theorem [7] holds for (|41|) if we require that b and o~i are continuous and bounded 
functions such that 

sup (\\b{t, -)\\ c e + \\o-i(t, -)\\ c i+e) < oo, i = l,...,k, 
te[o,T] b b 

and, moreover (as in Hypothesis [3|) we assume that o~(t, x) is uniformly non-degenerate, 
i.e., there exists the inverse of a(t, x) = a(t,x)o~*(t,x), for any t € [0, T], x € M d , and 

1 1 tz. 1 1 1 o = sup ||a~ 1 (t, x)|| < oo. (42) 

To prove Theorem [7] under these hypotheses, one can follow the proof of the analogous 
result proved in [10]. We only give a sketch of the argument. 

First note that |10[ Theorem 2] remains the same even with the previous non-constant 
a = (<7j) (indeed it is a special case of a result in [15]). Then [TUJ Lemma 4] is true with a 
in (I41j) by the following rescaling argument. Consider A > 1 and 

d t u\ + Lu\ — \u\ = f in [0, oo) x R d , 
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where L is the Kolmogorov operator associated to the SDE, i.e., 

L = -Tr[a(i, x)D 2 u{t, x)] + b(t, x) ■ Du(t, x) 

(here (cr(t, x)a*(t, x)) = a(t, x) and D and D 2 denote spatial derivatives). Define a function 
v on [0, oo) x M. d such that v(Xt, yXx) = u\(t, x), t > 0, x G M. d . It is easy to see that, for 
any s > 0, y G M d , 

M», ») + ^.(2, ^) »)] + ^(i • jt) ■ Dv{8 - v) - v{ '< v) = T f (i • Tx) ' 

Now the spatial Holder seminorms of (s,y) \— > a(j, and (s,y) h- > 6(j, are clearly 
independent on A > 1 and on s > 0. By [181 Theorem 2.4], we deduce in particular, for 
any A > 1, 

C 

sup||2>u(s,-)||o < T sup||/(s,-)||e, 

s>0 A s >o 

where C is independent of A. It follows the assertion of [101 Lemma 6] since 

C 

sup \\Du x (t, OHo = v / Asup||D?;(s,-)||o < -7= sup ||/(s, -)||e- 

t>0 s>0 vA s>0 

The proof of [10l Theorem 5] (which deals with the stochastic flow) remains true even with 
<T in (|41|) by a straightforward modification. 

Remark 10 An analogous of Theorem [7] holds for (|4ip requiring that Hypotheses [TJ [2] and 
[3] are satisfied "uniformly in time" . 

One assumes that b and Oi are continuous functions defined on [0, T] x M. d , i = 1, . . . , k. 
Moreover, there exists 9 G (0, 1) such that 6(t, •) G C e (R d ;M d ), i G [0,T], and sup tg[0 T ] \\b(t, •)||c e (R d ,»' i ) < 
00. In addition, ^(t, •) G C fe 3 (R d , IR 6 '), i G [0,T], 

SUp \\<Ji{t,-)\\ C 3 {R d Rd) < 00, 

te[o,T] 6 

z = 1, . . . ,k, and one requires that condition (|42|) holds. Theorem [7] under these assump- 
tions may be established by adapting the (time-independent) proof given in the present 
paper. However, the complete argument, even if it does not present special difficulties, is 
considerably longer (for instance, one has to prove the analogous of the Bismut-Elworthy-Li 
formula (|12p in the time-dependent case). 

We close the section by an application of the stochastic flow. We obtain a Bismut- 
Elworthy-Li type formula for the derivative of the diffusion semigroup (Pt) associated to 
([1]) (compare with [3] and [6]). It seems the first time that such formula is given for diffusion 
semigroups associated to SDEs with coefficients which are not locally Lipschitz. 
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Theorem 11 Let f : M. d — > R be uniformly continuous and bounded. For any x, h £ R d , 
we have (cf. O ) 

1 f 4 

D h P t f(x) = -E[f{<j> t (x)) / ({^a-^M^DhMxldWu)], t > 0, x £ R d , 
£ jo 

where (DP t f(x),h) = DhP t f(x) and Dcj) u (x) solves (f39j) anft s = (we sei <f> u {x) = 

Proof. We prove the formula when / £ C%°(M. d ). Indeed, then, by a straightforward 
uniform approximation of /, one can obtain the formula in the general case. 

Let i? : R d -> R be a smooth test function such that < &(x) < 1, x £ R d , = 
0(-s), J Rd i?(a?)da; = 1, supp (0) C B(0, 2), = 1 when x £ 5(0, 1). For any n > 1, let 
#n 0*0 = n dr d{nx). Define b n = b* # n . 

We have that b n is a C°° and Lipschitz vector field such that b — b n £ C^(R d ;R d ) and 
|| 6— 6 n |lc e tends to as n — > oo. Let (</>") be the associated flow of smooth diffeomorphisms 
which solves the SDE involving b n and let (P™) be the corresponding diffusion semigroup. 
The Bismut-Elworthy-Li formula for (P t n ) is given by 

D h P?f(x) = ±E[f(<t>?(x)) f\{a*a~ l ){ru{x))D h ci> n u {x),dW u )l t > 0, x £ R d , n £ N. 

Note that DhP™f{x) = E[(Df((pt(x)),Dh<fit(x))]. Passing to the limit as n — > oo, using 
the estimates (|29|) and ([3T]) . we get 

D h P t f(x)=E[(Df(Mx)),D h Mx))] = -nf(Mx)) / (a" 1 ^^))^^^),^)], 

1 Jo 

for any t > 0, x £ M d . ■ 
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